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Abstract 

A common experimental setup in cavity quantum electrodynamics (QED) consists of a single two- 
level atom interacting with a single mode of the electromagnetic field inside an optical cavity. The 
cavity is externally driven and the output is continuously monitored via homodyne measurements. 
We derive formulas for the optimal rates at which these measurements provide information about (i) 
the quantum state of the system composed of atom and electromagnetic field, and (ii) the coupling 
strength between atom and field. We find that the two information rates are anticorrelated. 
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I. INTRODUCTION 



In this paper, we consider a system consisting of a single two-level atom, located inside 
an externally driven optical cavity. The atom interacts with a single light mode inside 
the cavity, which is coupled to the environment via a partially transparent mirror. The 
output field from the cavity is monitored using the homodyne detection scheme, in which 
the cavity output is added to a reference field at a beam-splitter and then analyzed by two 
photodetectors (see Fig. HJ). 

nn 

In two recent ground-breaking experiments [H, 12], a similar setup was used to observe 
the trajectory of the atom inside the cavity. In these experiments, the atomic position was 
inferred from the strength of the atom-cavity coupling, which can be estimated directly 
from the photocurrents |3(. To characterize the performance of this atom-cavity microscope, 
Kimble introduced a quantity called optical information 0,0], which measures the rate at 
which the measurement provides information about the system. In Refs. , however, 
no formal definition of optical information is given, and only a heuristic derivation of its 
value is provided. Recently, a number of related definitions of information gain have been 
proposed and investigated numerically for a simpler quantum-optical system p. 

The question of how much information about a monitored system is provided by a contin- 
uous measurement is interesting in its own right. In this paper, we consider two, as it turns 
out, complementary cases. In Sec. IIII Al we calculate the optimal rate at which a homodyne 
measurement provides information about the quantum state of the system composed of the 
electromagnetic field and the internal state of the atom. In Sec. 1111 Bl we calculate the opti- 
mal rate at which the measurement gives information about the coupling strength between 
the atom and the intra-cavity field. In Sec. |H1 we introduce our mathematical model and 
main assumptions, and Sec. |IV| concludes the paper with a short discussion. 

II. MODEL 

The evolution of the state, p, of an open quantum system subject to a continuous mea- 
surement can often be described by a stochastic master equation of the form 

dp = C(p)dt + M(p)dW , (1) 



2 



which is understood in the sense of the Ito stochastic differential calculus jlo| . Any particular 
measurement record is represented by some realization of the stochastic process W(t). The 
superoperator £ is linear and defines the "unconditional" evolution, i.e., the evolution in 
the absence of measurements. By contrast, the superoperator M. is nonlinear and accounts 
for the effects of the measurement. 

The cavity-QED system we are considering here is illustrated in Fig. ^ A single two- 
level atom, with ground state |g) and excited state |e) is located inside a high-finesse optical 
cavity, which is driven by an external laser field of strength E. The atom interacts with a 
single light mode of the cavity. We denote by g the strength of the atom-cavity coupling, 
and by k the cavity field decay rate. In this paper, we assume that the atom, the cavity and 
the driving laser are all resonant. 

In the absence of measurements, M. = in Eq. and the joint density operator of 
atom and intra-cavity field, p, obeys the master equation 



p = £{p) 



(2) 



where 



C(p) = [E(a* — a) + g(a f a — era), p] + K(2apa) — a)ap — pa) a) . 



(3) 



Here a is the annihilation operator for the cavity field mode, and a = |g)(e|. In the following 
we focus our attention on the important case of the strong driving regime {E/g ^> 1). In 



this regime the system approaches a steady state of the form 

1 
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-[\a;+)(a;+\ + \a*;-)(a*;-\) 



Here \a; +) and \a*; — ) are two orthogonal quantum states defined by 

!";+> = -^|a)(|g) +*») , 

!«*;-> = ^K)(|g) -*>)), 

where \a) is the coherent field state with amplitude 



E 

a = — 

K 



\2EJ 2E\ \2EJ 



(4) 



(5) 



(6) 



Let us now assume that the field escaping from our system undergoes a continuous mea- 
surement via the standard homodyne measurement scheme In this scheme we have 




FIG. 1: Homodyne or heterodyne measurements in cavity QED. See the text for an explanation of 
the parameters. 

only one free complex parameter: the reference field /3 that is added to the cavity output at 
a beam splitter prior to the detection by two detectors D\ and D 2 (see Fig. GJ. If the mea- 
surement record consists of the scaled difference photocurrent dq/dt = (7 2 — h)/\P\, where 
l\ and I2 are the photocurrents detected by D\ and D 2 , respectively, then the measurement 
operator M. in the stochastic master equation (JTJ) is given by 



a 



M (p) = (e^ap + e i(t> pa ] - tr [p{e~^a + e^a f )] p) , (7) 

where rj is the efficiency of the photodetection and <p = arg f3 is the phase of the reference field. 
The realization of the Wiener process W(t) in Eg. O is connected to the experimentally 
observed difference photocurrent via the relation 

m 

dq = 2kt] tr [p(e*V + e -< *a)] dt + dW . (8) 

In the special case p = p" s this equation becomes 



dq = 4kt] Re (a) cos (j) dt + dW . (9) 
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III. INFERENCE 



A. Information about the quantum state 

The amount of information provided by the measurement about the quantum state, p, is 
quantified by the reduction of the von Neumann entropy, 

H(p) = -tr(plnp). (10) 

In this subsection we calculate the rate at which this entropy changes as a result of mea- 
surements. The average rate of entropy change in the presence of continuous observations 
is given by 

where the average (•) is taken over all possible measurement outcomes observed over the 
time At, and 

5 = £{p)At + M(p)AW . (12) 
We now assume that the system is initially in a steady state, p ss , defined by the relation 

£(Pss) = . (13) 

For simplicity, we furthermore assume that 

[M( Pss ), Pss ] = . (14) 

This last assumption is satisfied by Eqs. (@J) and ((7j). At the end of this section we will 
outline a more general framework which is free from these assumptions. Since the stochastic 
master equation ((TJ) preserves the positivity of the density operator, p, the support Ai(p ss ) 
must lie in the support of p ss , supp(A^(p ss )) ^ supp(p ss ). Now consider a basis in which p ss 
and A4(p ss ) are both diagonal. In this basis, let a\, . . . , a m be the nonzero diagonal elements 
of p ss and let bi, . . . ,b m be the corresponding diagonal elements of Ai(p ss ). 
Substituting p = p ss into Eq. (fTTj) . we obtain 

(HM) = ^( Pss + M( P ^AW)-H( Pss ^ 

lim / ~ EZLifofc + h AW) \n(a k + b k AW) + Eg=i ^ lna fc \ 
A™0\ At I 
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E 

fc=i 



2a k 



(15) 



where we have expanded the logarithm to second order in AW and used the relation 
((AW) 2 ) = At. 

For any state p° of the form (J3J we have 



M(pg S ) = sin0 \/2kt] Im(a) +)(a; +| — |a*; —)(«*; 



(16) 



Since p" s and .M(p" s ) commute, we can apply Eq. (jTKjl . For the particular value for a given 
by Eq. ©, the diagonal elements of p" s and M.(p^ are 



ai, 2 = l/2, and 6i j2 = ±# sin J ^- (l - (^ 



(17) 



Finally, we obtain the main result of this subsection: the rate of information gain about the 
quantum state of the system is 



(18) 



Here the minus sign indicates that the information gain corresponds to the reduction of 
uncertainty about the system state which is measured by the von Neumann entropy. 

We now finish this subsection by developing some general formalism that makes no as- 
sumptions regarding the superoperators £ and M.. Equation (JTUJ) can be expanded in the 
form 



H(p) = t r[^^Lp(p-ir 

| Z — / n 



tr 



n=l 

oo 

E 

n=l 
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n 



-(p n+1 + p n ) 



(19) 



where p = p — 1. Equation (jll|) then becomes 



•l) n 



^o°)> = j&E ^at <tr[G5 + 6)n+1 ~ ~ pn+l] + tl[Cp + 6)n ~ pn]) 

n=l 

Keeping the terms to the second order in 5, we have 

n— 1 n—2 n—2—p 

(p + S) n = p n + J2p l 6 P rM + P" 6 p n ~ 2 ~ ip+Q) 5 p q + 0(5 3 

1=0 p=0 9=0 



(20) 



(21) 
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Using the cyclic property of the trace, 

n—2 n—2—p 

tr[(p + 5) n - p n ] = tr Inp 11 - 1 5 + ^ P^ 9 6 p n " Mp+<7) <5 + 0{5 3 

n-2 



tr 



np 



5 + J2( s + 1 )P S 5 p n ~ 2 ~ s $ + 0(5 3 ) . (22) 



s=0 



Combining Eqs. (J2UJ) and ()22jl. we obtain 

oo /_-|\n , 

= ^AT tr ( (n + l) ~ pn5+ S (s + V? 6 ?^ 1 " 6 

n=l ' s=0 

n-2 

+ np n - 1 5+ + l)p s Sp n - 2 ~ s 5 ) . (23) 
Substituting Eq. (JT^J) and using (AW) = 0, ((AW) 2 ) = At we finally arrive at 
(H(p)) = ^-Ur( [(n + l)p n + np n - l ]C(p) 

n=l 

n-1 

+^( s +l)p s >f(p)^- 1 -^(p) 

s=0 

n-2 

+ J> + l)f?M(p)/r-*-M(p) ) . (24) 

s=0 

This formula is valid without any restrictions on the superoperators L and A4. It is straight- 
forward to show that, under the previous assumptions p = p ss and [A4(p ss ), p ss ] = 0, Eq. (|2*4*j) 
reduces to Eq. (|T5j) as required. 

B. Information about the atom-cavity coupling 

Since the quantized degrees of freedom described by the density operator p of the pre- 
vious sections do not include the atomic position, a different approach is needed to obtain 
information about the atom's trajectory inside the cavity. The key is the fact that the atom- 
cavity coupling g depends on the position of the atom in a known way (see, e.g., Ref. 0|). 
In this subsection, we will therefore focus on calculating the rate at which the measurement 
provides information about the parameter g. 

Let q be a photocharge obtained by integrating the difference photocurrent, given by 
Eq. (jHJ), over a small time interval At. Introduce a small parameter e such that nAt ~ e 3//2 
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and \(3\ ~ e -1 . The existence of such e is assumed in the standard derivation of Eqs. (J7J) 
and (JHJ) j^, |^ . Equation (JHJ) is derived in the limit of small e when q can be treated as a 
Gaussian random variable, 

fi ('M S ^(- ! Sf), (») 

with mean 

(g) = 2/^ tr[p(e^a f + e~^o)] At + 0(e 2 ) (26) 

and variance 

v 2 = 2n V At + 0(e 2 ) . (27) 

So far we have treated the atom-cavity coupling g as a known parameter. This means 
that the above equations give us the conditional probability density P{q\g) of registering 
measurement result q given that the atom-cavity coupling is equal to g, 

P(q\g) = G(q,(q),v 2 ). (28) 

The conditional probability of the atom-cavity coupling to be equal to g given a particular 
measurement result q can be derived from the Bayes rule, 

P(q\q) ~ (29) 

where P(g) is a probability distribution that characterizes our knowledge of g prior to 
obtaining q. 

To quantify the information gain, we use the entropy S of a continuous probability dis- 
tribution f(x) defined by 



S[f(x)] = -J f(x)\nf(x)dx. (30) 
The average rate at which observation of q gives us information about g is then given by 

«»-i;s„(^) = Ii n ;„i/(- A5 ) p («^' < 31 > 

where 

AS = S[P(g\q)] - S[P(g)\ . (32) 

On the relatively slow time scales of the atomic motion one can assume to a very good 
approximation that the quantized degrees of freedom described by the density matrix p are 
in the steady state given by Eq. (jlj). In this case Eq. (j2T?j) gives 

(q-(q)) 2 = q 2 -4qK V tr[p{e i(i, al + e-^a)} At + O (qe 2 ) + O (e 3 ) 
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= q 2 - 8qEr] cos (l - (^) 2 ) At + 0(qe 2 ) + 0(e 3 ) . (33) 

We can now derive a simple formula for R g if we make the further convenient assumption 
that the prior distribution P(g) is a Gaussian with variance v 2 , 

P(g) = G(g,(g),v 2 ). (34) 

The posterior probability density P{g\q) will then also be a Gaussian, 

P(g\q) = G(g,m,v 2 1 ) , (35) 

where the variance v\ is given by 

«? = zt 1 VlKE 2 1 ■ (36) 
kE + gt> o cos 

This expression, where we neglected the 0(qe 2 ) and 0(e 3 ) terms in Eq. (|33jh is valid for small 
g. The value of the mean m is irrelevant for the current argument. By a direct calculation, 
we obtain the entropy change 

AS = -- ln(v /vf) = -- In (1 + q-^v ) ~ -<l7^ v o ■ ( 37 ) 

Due to the limit At — > 0, the Gaussian P(q , |(7) in Eq. ()3H) is strongly peaked at small values 
of q. Hence we can substitute the above estimation of AS into Eq. (f3*Tf . We thus obtain 
our final result, the average information rate 

^Ti'-fslV*' (38) 

IV. DISCUSSION 

Comparing the expressions Eq. (|18|) for the information rate Rq and Eq. (|38|) for the 
information rate R g , one notices their very similar structure. The rate Rq, referring to 
information about the quantum state of the atom-cavity system, is proportional to sin 2 (f>, 
whereas the rate R g , referring to information about the parameter g, is proportional to 
cos 2 4>. In both formulas, (f) is the phase of the reference field (3. The proportionality factor 
in R g is obtained from the proportionality factor in Rq by substituting for g 2 twice the prior 
variance, 2v$, of the random variable g. 
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This means that there is a trade-off between the two information rates: the more we learn 
about the coupling parameter (and hence the atomic position), the less we can learn about 
the atom-cavity quantum state and vice versa. This is a manifestation of the uncertainty 
principle applied to the conjugate field quadratures X = a + a* and Y = i(a — a'). Let us 
elaborate on this point by considering our starting equations, i.e. Eq. (J8J), which gives our 
measurement record, and Eqs. (JTJ), (J3J), and (JTJ) for the corresponding change of the system 
density matrix. As before, we assume the steady state p ss as the initial condition for these 
equations. The qualitative conclusions of this discussion do not depend on any specific prior 
distribution for g. 

We notice that the average photocurrent is proportional to the expectation value of the 
expression e^a* +e~ lr ^a. In this sense the cases = and = tt/2 correspond to measuring 
the expectation values of X and Y quadratures respectively. In the case of = the 
nonlinear contribution from the measurement, given by Eq. (JTJ), becomes proportional to p" s . 
From Eq. (JJJ we therefore see that the system density matrix remains equal to p" s at all times, 
independent of the measurement record. This means that for = the measurement record 
contains no information about the quantum state of the system, in agreement with (|18j) . 
At the same time the setting = maximizes the average photocurrent, which carries 
information about g [see Eq. ©]. The second stochastic term in (JHJ) does not depend on 
g and so, from the point of view of learning about g, this stochastic term is nothing but 
noise superimposed on the average photocurrent. We therefore see that maximizing the 
average photocurrent increases the signal-to-noise ratio for measuring g and thus = 
gives us the maximum information about g, in agreement with Eq. (JHBJ)- In summary, the 
case = is ideal for measuring g and it gives no information about the internal state. In 
the case = it/2 we have the opposite situation: we obtain maximum information about the 
atom-cavity state and no information about g, because the measured photocurrent becomes 
independent of g. For a numerical comparison of the two information rates for a simpler 
quantum-optical system and a wide range of detection schemes, see Ref. jfj|. 

At first sight, learning about the quantum state of the system seems to be a very different 
problem form learning about the atomic position. It so happens that within our framework 
these two tasks are best accomplished by monitoring conjugate observables, namely the X 
and Y field quadratures. With a straightforward modification our calculations can be also ap- 
plied to the heterodyne detection scheme. For heterodyne detection, both field quadratures 
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are monitored simultaneously, at the cost of reducing the respective maximal information 
rates. 
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